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Abstract 



We study the vortex state in a magnetic field parallel to the c axis in the 
framework of the extended Ginzburg Landau equation. We find the vortex 
acquires a fourfold modulation proportional to cos (4</>) where eft is the angle 
r makes with the a-axis. This term gives rise to an attractive interaction 
between two vortices when they are aligned parallel to (1,1,0) or (1,-1,0). 
We predict the first order vortex lattice transition at B = H cr ~ K~ 1 H C 2(t) 
from triangular into the square lattice tilted by 45° from the a axis. This 
gives the critical field H cr a few Tesla for YBCO and Bi2212 monocrystals at 
low temperatures (T < 10K). 
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The o?-wave superconductivity in the hole-doped high T c cuprate superconductors now 
appears to be universally established fl],|2[]. Then in a magnetic field parallel to the c axis we 
expect a square vortex lattice tilted 45° from the a axis at least in the vicinity of B = H c2 (T) 
and below t(= T/T c ) < 0.81 Indeed such vortex lattices though elongated in the a di- 

rection have been seen by small angle neutron scattering (SANS) H and scanning tunnelling 
microscope (STM) |J both in twinned YBCO monocrystals at low temperature (T ~ 4K) 
and in a small magnetic field (B ~ a few Tesla). The object of this letter is to determine 
theoretically the corresponding vortex lattice transition at a small magnetic field which is 
experimentally relevant. Also in order this vortex lattice transition is to be observable ex- 
perimentally, it should take place in a magnetic field smaller than the one where vortex 
lattice melting is established Further we shall show that the transition is exthermal 

and of the first order. Earlier the related transition has been analyzed in terms of the 
generalized London equation, which may be obtained from the phenomenological Ginzburg 
Landau equation for c?-wave superconductivity || or from the rhombic term arising from the 
normal state properties [ 10| . However these authors concentrated on the modification of the 
magnetic interaction energy due to the quartic term in the kinetic energy. We believe, on 
the contrary, that the important modification lies in the vortex core structure [0]]. Never- 
theless these authors obtained the vortex lattice transition from a triangular lattice in small 
magnetic field to a rhombic lattice in a higher field. The result is similar to ours except for 
two crucial differences. First we find that the critical field H cr where this transition takes 
place behaves like K^H&it) where k is the Ginzburg Landau parameter and H C 2(t) the 
upper critical field. This means H cr ~ 10~ 2 H C 2(t) in high T c cuprates. Second the transition 
is of the first order with a sharp negative divergence in the magnetization. 

Within the weak-coupling model for <i-wave superconductor |12|] the extended Ginzburg- 
Landau equation is given by |TT] 
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where t = T/T c . Here d x and d y are the gauge invariant differential operators. Eq.(^) has 
been obtained previously by Enomoto et al. fL3| , though we suppressed all the other terms 
of the order of (— lnt), which is not necessary for our purpose. Due to the quartic terms in 
Eq.(|l]) the single vortex solution is now given by [|ll|] 
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where g(r) ~ tanhr and |/(r)| and the phase of f(r) are shown in Fig.l and Fig. 2. For 
clarity we have chosen e = 1. The ridges in |/(r)| running along four diagonals are clearly 
seen, which are crucial for the fourfold core interaction energy. Also the phase distortion 
propagates far away from the vortex center. For r 3> 1, a(r) and (3{r) are obtained as [|TT 
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which tell us the explicit form of the asymptotics. This result is consistent with Enomoto et 
al. ||13|| . In the presence of the fourfold distortion around single vortex, this term modifies 
drastically interaction between two vortices. Now let us consider the vortex lattice. Within 
the dimensionless unit the free energy of the vortex lattice is given by; 
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where b(r) is the local induction. Making use of the usual approximation for A(r) = 
^I\if{ r ~ r i) where is the position of the 2-th vortex. We are here interested in the 
region 1 <C d <C k (or £(T) <C d <C A(T) in the natural unit), where d is the distance 
between 2 neighboring vortices. Then neglecting terms exponentially small we obtain 



n 



£4 

A — ai^n^ — elOai^n^ cos 49, 



Lm r l,m 



Lm 



2tt 



K 



Lm 



n,rr 
X 



where A is the area, a% = ^ L (ln2 + |) — 6.354, = B/<p the vortex density per 
unit area and r^ m = 2ld(cos 9 , sin 9) + 2md(cos9, — sin 9) or (21 + l)d(cos9,sin9) + (2m + 
l)d(cos9, — sin 9). Here we limit ourselves to the vortex lattice with isoceles unit cell. Also 
as to the magnetic interaction we consider here only the isotropic term k~ 2 K (j^) since 
the anisotropic correction is much smaller than the vortex core interaction at least when 
k 1. The last term in Eq.( 
we transform Eq.(RI) as 
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is handled using the Poisson transform as in [fL4|| . Finally 
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By minimizing f(9) in 9, we find the apex angle 9 B 
gradually and jumps to 7r/4 at B = H cr where 



shown in Fig. 3. 9 min increases first 



H cr = 0.524(-lnt) 



-1/2, 
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Earlier a similar 9-B curve is obtained within generalized London equation P/TCfl . But here 



the arise of 9 to 7r/4 is much steeper. Indeed contrary to H[TIJ, we predict that the transition 
is of the first order with a small dip in — M [p"T . 



Unfortunately no experimental data on 9 m i n is available in the hole doped high T c cuprates 
to our knowledge. However a very similar jump in 9 m i n as function of B as described here 
has been seen by SANS from ErNi2B2C at T = 3.1 jl5fl . Although the superconductivity 
of borocarbides is believed to be of s-wave, both the existence of the rhombic vortex lattice 



in a magnetic field B parallel to the c axis seen by SANS [p~5|JT6|1 in ErNi 2 B 2 C, YNi 2 B 2 C 
and LuNi 2 B 2 C and STM from YNi 2 B 2 C [17] suggest strongly that the underlying supercon- 
dudtivity is of <i-wave. More recently the fourfold anisotropy in the uppercritical field in a 
magnetic field within the a-b plane seen in LuNi 2 B 2 C is again consistent with rf-wave 
superconductivity |19| , |20[ . 



To sum up we find that the vortex lattice is triangular as in a classical s-wave supercon- 
ductor in a small magnetic field. As B increases the vortex lattice transforms first gradually 
and then rapidly into a square lattice at B = H cr (t). Also the final jump is extremely rapid, 
which is seen from a sharp dip in the magnetization. 

Although the theoretical expression for H cr (t) is obtained within the weak-coupling model 
and for t not too small (say t > 0.5), from the temperature dependence of C{t) in 0, we 
expect the result should be valid semiquantitatively even for t <C 1, if we replace — hit by 
1 — t. This allows us not only to estimate H cr (0) in the optimally doped YBCO and Bi2212 
crystals but also to construct the phase diagram for the rhombic vortex lattice. As readily 
seen from Fig. 4 the rhombic lattice occupies the major part of the T-B phase space of the 
vortex state unless disturbed by the vortex lattice melting. 

Though -ff c2 (0) in YBCO or Bi2212 is not known for sure, we may assume if c2 (0) ~ 120T 
and 300T for YBCO and Bi2212 respectively. Then making use of the fact k ~ 100 for these 
systems, we expect that the vortex lattice transition takes place at B ~ IT and 3T for 
YBCO and Bi2212 respectively. The former is consistent with SANS and STM experiment 
in YBCO monocrystals mentioned in the beginning. For Bi2212, on the other hand, the 
melting transition appears to take place around B ~ 300 Gaufi for optimally, doped Bi2212 
||. If it is really the case, it will be difficult to see the square lattice discussed here in Bi2212 
unfortunately. On the other hand if this transition exists, it will be not difficult to identify 
by SANS, STM and even by thermodynamic measurement. 
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Figure Captions 



Fig.l: |/(r)| in the x-y plane, we used e = 1 for the clarity of the figure, while e ~ 0.1 in 
the real situation. 

Fig.2: $(r) = phase of /(r) — is shown in the x-y plane for e = 1. 

Fig.3: Appex angle 26> min as a function of B/H cr where 29 min = 90° and 120° correspond 
to the square lattice and the triangular lattice with hexagonal symmetry, respectively. 

Fig. 4: The phase diagram with the rhombic vortex lattice is shown schematically. In 
reality the rhombic lattice becomes stable around B = 10~ 2 H c2 (t). Also "irr" means 
the irreducible line where the vortex lattice melts into the vortex liquid by the first 
order transition in the clean high T c cuprate superconductor. 
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Fig.l 




Fig.3 



